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ABSTRACT 

A straightforward analytic approximation is illustrated to 
Stromgren's solution for HII regions which enables us to explicitly 
exhibit the effects of internal dust on the ionization structure. 

We interpret recent far-infrared observations of HII regions in terms 
of true absorption by internal dust of a significant fraction of the 


Lyman continuum photons , 



The discovery of large far-infrared emission from HII regions- 
(Harper and Low 1973-) > significantly in excess of the apparent Lyman 
continuum output of the exciting stars inferred from measurements of 
the thermal radio continuum, on the assumption that the Lyman con- 
tinuum photons are absorbed by the gas has stimulated us to re-examine, 
the effects of internal dust on Stromgren's classical solution to the 
HII region problem. It 'is the purpose of this note to indicate a 
simple analytic approximation, which leads to results in reasonable 
agreement with earlier numerical work (Mathis 1973-) • Our result more- 
over both clarifies and adds considerable insight into the problem, 
and is useful in indicating how one may extend Mathis' results over 
a wide range of optical depths for dust absorption, with no additional 
numerical computation. Moreover, we apply our results to interpret 
the observed far-infrared emission from HII regions, and demonstrate 
the observations to be consistent with a moderate amount of absorption 
of Lyman continuum photons by internal dust. 

The standard equation of ionization equilibrium for a spherically 
symmetric HII region of uniform density n, containing a central 
exciting star emitting S ionizing quanta per second can be written 
in the form 
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where x is the fractional degree of ionization, y = r/r g the 

dimensionless radial coordinate in units of the Stromgren radius, 
which is r , defined by 
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In these equations CH is the hydrogen recombination coeffi- 
cient summed over all levels excluding the ground level, and a is 
an appropriately averaged value of the photo-ionization cross-section 
for hydrogen. The optical depth t is defined by 
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where 
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cr is the effective absorption cross-section of the grains, and x 
g g 

is the number of dust grains per H atom. The dust and gas are assumed 
to be uniformly distributed in the ionized region. 

We make the approximation 1-x « 1, which should be well satisfied 
throughout the ionized region. It is then a straight-forward matter 
to eliminate T from equations (l) and (4), and obtain the following 
solution for x as a function of radial coordinate, 
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For « 1, equation (6) reduces to the anticipated limit 
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The neutral fraction 1-x is plotted in Figure 1 as a function 
of y for various values of - The parameter T g is set equal 

to 5000- F° r comparison, the nxunerical results of Mathis ( 1971 ) are 
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also shown (dashed line) for one of the cases he computed, namely 

T = 0.4. The agreement for the ionization structure is remarkably 
d 

good, despite the approximations of our analysis. The small disagree- 
ment at low values of y is due to scattering of photons by dust 
which we have ignored completely. 

Another useful quantity that we can compare with Mathis' numeri- 
cal work is the fraction of Lyman continuum photons absorbed by the 
gas. The analytic expression for this fraction as a function of y 
is found to be y J /(l-e ) . At the boundary of the HII region (y = y ) 
(1-x)t » 1, and t -* «=. Consequently, we obtain the extent of the 

ionized region by setting the left hand side of equation (7) equal to 
zero. The fraction f of stellar ionizing photons absorbed by gas 
in the entire ionized region is therefore given by 

f = y Q 3 = 3 T 0 3 e 0/ [ t 0 2 “ 2t 0 + 2(1-6 ’ (8 ) 

where n T^y^ is the optical depth for absorption by dust within 

the ionized region. 

In Figure 2 we plot the fraction f versus t . For comparison 

with our analytic curves, we show results computed numerically by 

Mathis. As is evident our solution is in perfect agreement with the 

computations of Mathis for t ^ 1.2. For larger amounts of dust 

o 

absorption, our solution tends to overestimate the fraction of ionizing 
radiation absorbed by the gas. 

Part of this discrepancy is due to our neglect of helium. Mathis 
(1971) shows that the He + and H + ionization zones tends to coincide 
in the presence of dust, because the H ionization zone moves inward 
with increasing as H-ionizing photons are absorbed by dust. One 

of us has also given an analytic expression for this effect (Petrosian 
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1972) . The principal effect of He on the hydrogen ionization struc- 
ture is to act as a converter of relatively hard He-ionizing photons 
to softer photons (each He recombination liberting ~ 0.9 H-ionizing 
photons), for which the cross-section for ionizing H is much higher 
(Hummer and Seaton 1964). Consequently the ionization of H is in- 
creased, which implies that at ^ 1 our neglect of He leads us 

to overestimate the neutral fraction, and therefore the fraction of 
ionizing photons absorbed by the gas . 

Part of the different between the analytic and numerical, solutions 
appears to be due also to our neglect of scattering by the dust, the 
effect of which is to provide a higher probability for absorption of 
photons by either gas or dust grains . We can account for the effect 

c 

of scattering by redefining all our optical depths as t(t) • = j an ds, 

*0 

where ds is an element of length measured along the actual scattering 
path of the photon. With this new definition of optical depth our 
results remain unchanged, but the points in Figure ( 2 ) from Mathis' 


numerical work move to the right approximately by a fact 


r O . 
ds 
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where r Q = r^y^ Is the Stromgren radius of the dusty nebula. This 
factor increases with increasing values of the dust albedo. 

In order to apply Figure (2) to observations of far-infrared 
emission from HII regions, we note that the expected infrared lumino- 
sity is 

-T ' 

L(IR) - L(Lya) + (l-f)<hv) Lyc S + L v< Lyc (l-e °) ( 9 ) 

where is avera 6 e energy of stellar Lyman continmim 

photons, L is the stellar luminosity below the Lyman continuum 

’ v< Lyc 

and Tq is the effective absorption optical depth of dust for photons 
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with hv < 13.6 ev. For a nebula which is optically thick to Lya 
photons (case B), the value of the Lya luminosity of the nebula 

(assuming that the Lyman Alpha can eventually escape from the nebula) , 

2 

is given by L(Lya) = — Sf, and can be inferred either 

from the observations of Balmer photons suitably corrected for 
reddening or from the radio continuum luminosity. If the spectrum 
of the ionizing stars and the observed infrared flux of HII regions 
are known, we can use equation ( 9 ) to determine the fraction f, and 
from Figure (2) the effective absorption optical depth of the dust. 
This procedure may allow one to empirically determined the far ultra- 
violet absorption properties of the dust. For example, in the case of 
the Orion nebula, the total extinction by dust near H(3 has been found 
to be I m .8 (Mathis I 97 O) ; this value is consistent with other indepen- 
dent determinations at optical wavelengths (cf. Munch and Persson 
I 971 ). On the other hand, we estimate that, in our earlier notation 
L(lR)/L(Lya) = 7 - 5 - Hence if we neglect the last term in equation ( 8 ) 

as might be the case for a hot central star or if t ' « t , and set 

00 

M Lye * l 5 ev, we find that f « 0.26. We then deduce from Figure 2 

that the true absorption by dust of Lyman continuum photons is about 

I m .7. It follows from these numbers that, for a ratio of total absorp 

tion cross-section in the Lyman continuum to that at of between 

2:1 and 3:1; the albedo at H must be about 0.6 ± 0.1. In principle, 

P 

one could now infer further physical characteristics of the dust from 
this result: however, the paucity of available data on optical 
constants at UV wavelengths does not appear to allow any definitive 
conclusion to be drawn at present. In the event that the last term 
in equation ( 8 ) is not negligible, similar calculations can be carried 
out if the spectrum of the central star is known. 
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In summary^ our analytic approximation (6) to the ionization 
structure of HII regions with a uniform distribution of dust should 
be useful in aiding future studies of line and continuum emission 
from HII regions. In particular, we have obtained a simple expression 
for the fraction of ionizing photons absorbed by the dust [equation 
(8)], a quantity which is consistent with the infrared observations 
provided that the dust absorption optical depths are typically between 
1.0 and 2.0 magnitudes. This approach offers an independent way of 
obtaining the total absorption optical depth to internal dust in HII 
regions, which can then be compared with the optical depth at H|3. 

If the stellar spectrum is known, this approach can be used for deter- 
mination of the albedo of the dust. 
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FIGURE CAPTIONS 


Figure 1 . Log^ (neutral fraction) plotted against radial distance 
in units of r .. Each curve is labeled by the corres- 
ponding value of . For comparison, the numerical 
computation by Mathis (1971) is shown for the case 
= 0.4 (dashed line) . 

Figure 2 . Fraction of ionizing photons absorbed by gas plotted 

against optical depth for true absorption by dust within 
the ionized region. The numerical results of Mathis 
(1971) are indicated by crosses. 
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